Dynamics of Phase Separation under Shear: A Soluble Model 
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The dynamics of phase separation for a binary fluid subjected to a uniform shear are solved 
exactly for a model in which the order parameter is generalized to a n-component vector and the 
large-n limit taken. Characteristic length scales in directions parallel and perpendicular to the flow 
increase as {t^ / Int)^^'^ and (t/lnt)^'''* respectively. The structure factor in the shear-flow plane 
exhibits two parallel ridges as observed in experiment. 



The dynamics of phase separation under shear has at- 
tracted considerable theoretical 0, experimental 
and simulational [^-^ attention in recent years. In the 
absence of shear, the dynamics of phase separation is now 
quite well understood Domains of the two equilib- 
rium phases are formed, and coarsen with time in a man- 
ner well-described by a dynamical scaling phenomenol- 
ogy with a single growing length scale L{t) which gen- 
erally grows as a power law in time, L{t) ~ t". The 
structure factor is spherically symmetric, with a maxi- 
mum at wavevector km ^ . For binary fluids, the 
exponent a takes different values depending on the dy- 
namical regime under study. In order of increasing time, 
there are 'diffusive' (a — 1/3), 'viscous hydrodynamic' 
(a = 1) and the 'inertial hydrodynamic' (a = 2/3) 
regimes [p|-pT|. The crossover between these regimes is 
determined by the fluid properties (viscosity, density). 
Here we will focus on the diffusive regime, in which hy- 
drodynamic effects can be neglected. In the absence of 
shear, phase separation is described by the Cahn-Hilliard 
equation for the order-parameter field 0(r, i), namely 

at^ = ~v2(v2(/) + 0-03). 



If a uniform shear flow is imposed in the x-direction, 
with shear direction y, the flow velocity is Vx — jy, 
where 7 is the shear rate. For an incompressible fluid, 
the Cahn-Hilliard equation acquires an advective term 
V • \7(f) = jydx4> on the left-hand side. Generalizing to an 
n-component vector order parameter, this gives 



dt4> + iydx4> = 



For a critical mixture quenched into the two-phase re- 
gion from the homogeneous phase, an appropriate initial 
condition is a Gaussian random field with zero mean and 
short-range correlations: {(j)i{r,0)(j)j{r' ,0)) — A6ijd{r — 
r')- 

In this Letter we present the first exact solution for a 
system phase-separating under shear, by solving (^ in 
the limit n — > cx). We obtain characteristic length scales 
Lx ^ - ityint)'/\ Ly ^ fc-i ~ {t/\nt)y\ and 
Lz ^ ^ (t/ \nty^^ . These are extracted from the 
structure factor, S'(k, t) which has four maxima, located 
at k = ±{kmx, —kmyi^kmz)- Thus the rotational sym- 
metry of the zero-shear limit is completely broken. In the 



kz — plane there are two maxima terminating two long 
parallel ridges as observed in experiment ||]. The struc- 
ture factor exhibits multiscaling [l^ , rather than simple 
scaling, but this is presumably an artifact of the large-rt 
limit, as in the zero-shear case At the end of the pa- 
per we conjecture how the results will be modified for the 
scalar order parameter appropriate to binary mixtures. 

In the limit n ^ 00, one can replace /n in (|l|) by 
its mean in the usual way, leading to a self-consistent 
linear equation. After Fourier transformation this reads 



dt dky 



-k2[k2-a(i)]0k 



(2) 



where is (any) one component of 0, and a(t) = 1— {4> ). 

The same equation was studied in a recent Letter by 
Corberi et al. |Q, where it was regarded as a 'self- 
consistent one-loop' approximation to the scalar version 
of (^. These authors, however, did not solve the equa- 
tion analytically, but rather integrated it numerically (in 
two space dimensions). As a result they were unable to 
access the asymptotic {t 00) behavior which is the 
main focus of the present work. In particular we do not 
find the oscillatory large-time behavior that Corberi et 
al. conjectured from their numerical results. 

Equation (||) can be solved via the change of vari- 
kz,t) (fca;, (7, fcz, r), where t — t and 
a — ky + jkxt- The left-hand side of (||) then becomes 
d<p\i/ dr, and the equation can be integrated directly to 
give, after transforming back to the original variables, 
(1) 0k(i) = 0k(O) exp /(k, t), where 



/(k, t) ^-{kl + klft ~ {{ky + Ikxtf - kl) 



-t- {kl + kl + {ky + -fkxtf) b{t) 
-2jkx{ky + jkxt)c{t) + j^kle{t), 



(3) 



with b{t) = J^dt'a{t'), c{t) = ^'^ dt' t' a{t') , and e{t) = 

The next step is to determine a{t) self-consistently, us- 
ing its definition: a{t) — \ — (0^) — 1 — ^(k, t), where 
5(k, i) = (0k(i)0-k(i)) is the structure factor. The 
momentum sum can be evaluated for large t using the 
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method of steepest descents. To simplify the analysis we 
make the following ansatz, which will be justified a poste- 
riori. Naive power counting applied to (||) suggests char- 
acteristic length scales Ly ^ t^^^, and ^ t^/^, 
where the dominant /ca;-dependence comes from the shear 
term, and that a{t) ^ t^^/"^. This suggests that only the 
kx terms multiplied by 7 in survive in the 'scaling' 
limit. In fact we will find that the naive power counting 
result is modified by logarithms, but the above conclu- 
sion still holds. Guided by the 7 = result, we make 
the ansatz a{t) ^ (Int/t)^/^ for i — > 00. Then, to leading 
logarithmic accuracy, h{t) - (ilnt)i/2, c(t) tb{t)/3, 
and e{t) t^h{t)/b. Inserting these results in (^), and 
making the change of variable 

ikx = u, ky^^v, k^ = ]J^w (4) 



gives 



(5) 



S'(k, t) = A exp y2— F{u, v, w) 

1 8 4 

F{u,v,w) = -— Uu + v)'' -w^) + —u^ + -uv^v^ 
bu ' 15 3 

2 

--w^(u2 + 3mw + 3w2) + w2-w*, (6) 
o 

where contributions to F which vanish as t — *■ cxd (at 
fixed u,v,w) have been dropped. 

The self-consistency equation for a{t) reads 



1 - a{t) 



d^k 
(2^ 



S{k,t) 



(27r)37i5/2 

f f 

X / dudv dw exp l2— F{u,v,w] 



(7) 



Since b'^/t grows like Int by assumption, the triple in- 
tegral can be evaluated by steepest descents for large t. 
The complete set of stationary points (u, v, w) is listed in 
table 1, together with their type (maximum, minimum, 
saddle point) and the corresponding values of F. The 
stationary points with w = are also stationary points 
of the two-dimensional (2D) theory, and their types in 
the (u,v) plane are listed separately. The w ^ struc- 
ture of S(k, t) is relevant for light scattering with the 
incident beam normal to the xy-plane. 

The maximum value of F, Fm — 23/60, occurs at four 
points in uuui-space (labeled / in Table 1), correspond- 
ing, via (W), to four points in momentum space. The in- 
tegral in (0) can now be evaluated by steepest descents. 
Since a{t) in (Q) vanishes (like (Int/t)^^'^) for i — > cx3 it 
can be dropped to give 



1 = const. 



A 



■ exp 



2K^ 
t 



whence, to leading logarithmic accuracy, 
finally justifying our original ansatz. 



(8) 



(9) 



Label 


Position 


Number 


F 


Value 


Type (3D) 


Type (2D) 


a 


(0,0,0) 


1 








Min 


Min 


b 


±(2/V3,0,0) 


2 


16/45 


.35556 


IS 


IMax 


c 


±(72- 1/V3,-1/V2,0) 


2 


(37- 12V6)/180 


.04225 


S2 


S 


d 


±(V2 + 1/V3,-1/V2,0) 


2 


(37+ 12V6)/180 


.36885 


SI 


Max 


e 


±(0,0,1/V2) 


2 


1/4 


.25 


SI 




f 


±(^,-l/\/3,±l/^) 


4 


23/60 


.38333 


Max 





Table 1. Stationary points of F{u,v,w): Max = maximum, Min — minimum, S = saddle point (2D), Sn = saddle 
point of type n (the matrix of second derivatives has n positive eigenvalues), IS = 'inflection saddle point' (one 
positive, one zero, one negative eigenvalue), IMax = 'inflection maximum' (one zero, one negative eigenvalue). 



From equation (Q) we can define characteristic length 
scales for the three directions: = 'y{t^/by^^ ^ 

7(iVlnt)i/4, and Ly ^ L, ^ {t/bf'^ - (t/lnt)i/4, 
by setting u = kxL^, v — kyLy, and w = kzLz. 
Finally the structure factor is given by 5'(k, t) — 
A exp[(26^/i)F(u, w, w)]. Using (H) this can be recast as 

5(k,t) - const. (lnF3)3/2-l//('i)/^'", (10) 

where Vs{t) = L^LyLz ^ -ft' /"^ / {Intf/^ is the 'scale 



volume' at time and q = (k^L^, kyLy, k^L^) is 
the scaled momentum. The corresponding result for 
zero shear (in space dimension d) is [p|jl^] S'(k, t) = 
const. (In L)i/2Ld'^'(9)^ where L = (St/dlnt)^/'^, q = kL, 
and 0(9) = 2g2 _ Equation (|l|) does not have a form 
consistent with conventional scaling, 5(k, t) = Vs{t)f{q), 
but rather a form of 'multiscaling', as in the zero-shear 
case, where the power of the scaling volume Vs depends 
on the scaling variables. As in the zero-shear case, how- 
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ever, we anticipate that simple scaling will be recovered 
asymptotically for any finite n, and that the Int terms 
will disappear from the length scales . 
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FIG. 1. Contour plot of F{u,v,0) showing approximate lo- 
cations of the stationary points a, b, c, d from Table 1. Contour 
lines for F < —0.1 are not shown. 

In measurements of the structure factor by small- 
angle light scattering, the scattering vector k is per- 
pendicular to the beam direction. Two common exper- 
imental arrangements |^ are with the beam in the z- 
direction, i.e. perpendicular to both the shear and the 
flow, or in the y (i.e. shear) direction. The former 
corresponds to fc^ = 0, the latter to ky — 0. For 
kz — 0, S{kx,ky,0,t) is determined by F{u,v,0). In 
fact, from (0), ln5(k,i) = [F{u,v,0)/ F^]\nVs (plus 
k-independent terms), so F{u, v, 0) is essentially the log- 
arithm of the structure factor divided by Int. Figure 
1 is a contour plot of F(u,v,0), showing the stationary 
points a, 6, c, d listed in Table 1. The two global maxima, 
labeled d, are connected by long, almost straight, ridges 
to the two 'inflection maxima', labeled b. Note that the 
maxima, at a 'height' of 0.36885, are little higher than 
the inflection maxima, at 0.35556, so the ridges are al- 
most level. This ridge structure is strikingly similar to 
what is observed in experiments Note that the rela- 
tions InS" cx Flnt, and k^ — u/L^ etc., mean that the 
ridges in S become higher, narrower, and closer together 
as a function of k with increasing time. We can calculate 
the angle, 0, of the ridges to the w-axis (shear direction) 
from the slope of the line joining points b and d , giv- 
ing tan0 = 2(1 — l/\/6) (the same result is obtained 
from the line joining the saddle points, c). Using (^, the 
corresponding angle, 9, in the {kx,ky) plane is given by 



tanO = (tan0)/7<, i.e. the ridges tend to align closer to 
the shear direction as t increases ||]. 

For the beam in the y direction, the scattering inten- 
sity, S{kx,0,ky,t), is determined by F{u,0,w). Figure 
2 is a contour plot of this function with the positions 
of the stationary points {u, w) indicated. These are a 
minimum (labeled a) with F = 0, at (0, 0), a pair of sad- 
dle points (5), with F — 16/45, at (±2/-\/3, 0), a second 
pair of saddle points (c), with F — 1/4, at (0, ±l/\/2), 
and four maxima (d), with F = 29/80, at (±3, ±1)/2V2. 
Currently available experimental data cannot resolve this 
structure. Instead a broadly elliptical scattering pattern 
is seen in the {kx, kz) plane with the major axis along 
the kz direction, and the eccentricity increasing in time, 
as expected from the different growth rates (by a factor 
t) of Lx and Lz- 
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FIG. 2. Contour plot of F{u, 0, w) showing approximate lo- 
cations of the stationary points a,b,c,d (see text). Contour 
lines for F < —0.1 are not shown. 

Finally, if the beam is parallel to the flow {kx = 0), 
the shear term drops out of (H). The scattering intensity 
then has full circular symmetry in the (ky, kz) plane. 

From the asymptotic solution, the expressions derived 
by Onuki for the shear-induced contributions to the 
viscosity. At] = — (I/7) /[d^fc/(27r)^] k^ky S{k.,t), and to 
the normal stress differences AA^i = /[d'^fc/(27r)'^] [ky — 
kl) S'(k, t), and AN2 = /[d3fc/(27r)3] (fc2-fc2) 5(k^ t) may 
be easily evaluated. Since 5(k, t) is a sharply peaked 
function, the factors involving kx, ky, kz in the integrands 
can be replaced by their values at the peaks (/ in Table 
1). Using also /[(i^fc/(27r)^] ^(k, i) = 1 gives, asymptot- 
ically, Ar/ = (V7^i^) (lnt/t3)i/2^-2^ ^ (5/3^) ^ 
(Xnt/tf^, and AiVa = (&/6t) = AiVi/2. 
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We now compare our analytical results with the nu- 
merical solution of (H) by Corberi et al. in two di- 
mensions (the shear-flow, i.e. xy, plane). Repeating our 
analysis for d = 2 is straightforward. The dominant sta- 
tionary points are the two maxima labeled d in Table 
1 and Figure 1. A two-dimensional steepest descent in- 
tegration over the scaling variables u and v yields the 
self-consistency condition 

1. const. -exp(^^j , (11) 

instead of (|), where (from Table 1) F^^ = (37 + 
12v^)/180, giving b{t) = (3tlni/4i^™^)^/2 (to leading 
logarithmic accuracy) instead of (|9|) . The structure factor 

is given by 5(k, i) = const. (In A^) , where 

As = LxLy is the 'scale area', — ^{t^ /bY^"^ ~ 
7(^Vln^)'/^ Ly = [t/hf^ ^ {t/\ntY'\ q - 
{kxLx,kyLy) is the scaled momentum, and F^^{q) 
means F{u,v,0), with u — kxLx and v — kyLy. In 
other words, the 2D structure factor has an identical 
shape (shown in Figure 1) to the 3D structure factor with 
kz — 0. The main difference is that the numerical values 
of Frn (0.38333 . . .) and F^^ (0.35555 . . .) are (shghtly) 
different. 

In their numerical results, Corberi et al. also find two 
ridge-like structures, but the ridges are terminated by 
two peaks whose relative heights oscillate in time, such 
that first one peak, then the other is the higher. They 
further speculate that these oscillations persist to late 
times and 'characterize the steady state'. Since no oscil- 
latory behaviour is found in the asymptotics of our ana- 
lytical solution, we believe that the observed oscillations 
are slowly-decaying preasymptotic transients. 

We conclude with some conjectures about the phys- 
ically realistic case (for binary fluids in the diffusive 
regime) of a scalar order parameter. These are informed 
by our exact solution for n = oo, and by the way the 
n = oo solution is known to be modified for scalar fields 
in the unsheared case [ p^ . First we expect that, for any 
finite n, the structure factor will exhibit asymptotic scal- 
ing of the form S'(k, i) = Vsg{q), with Vs = Ili=iLi 
{i — x,y,z) and qi — kiLi, instead of the multiscaling 
form (p^. As in the n = oo case, we expect the growth of 
the characteristic scales for directions normal to the flow 
to obey the same power laws as in the unsheared case, 
i.e. for scalar fields, Ly ^ t^^^. The growth in the 

a;-direction can then be deduced from the assumed scal- 
ing form for the structure factor: if we multiply the two 
terms on the left-hand side of (||) by (/)_k(i)j and average, 
the result ^ itLy follows immediately if we insert the 
scaling form for ^(k, t) and assume both terms are of the 
same order in the scaling limit. This leads to the predic- 
tion Lx ~ 7<*^'^. The shear-induced viscosity and normal 



stresses then scale as A77 ~ (pjLxLy) ^ ~ 7 "^t and 

In the viscous hydrodynamic regime, where without 
shear L{t) t, the same heuristics suggest Ly ^ ^ t, 
Lx ~ (with corresponding modifications to A?7 and 
AiVi.2). The predictions for the length scales are consis- 
tent with data on polymer blends , though it has been 
suggested ||l|,^ that a stationary state eventually devel- 
ops for 7t 3> 1 due to a competition between stretching 
and breaking of domains. 

As far as the shape of the structure factor is concerned, 
we expect the scaling function ^(q) to be described, in 
broad terms, by ln(7(q) ~ F{qx,qy,qz), where F{u,v,w) 
is given by (^), i.e. rather like the large-n theory but 
without the InV^s ^ Inf factor multiplying F. In par- 
ticular, a contour plot of \ng{qx,qy,Q) should look very 
similar to Figure 1. (However, the behaviour near g = 
would be modified by the requirement that g(Qi) = 0, im- 
posed by the conservation of the order parameter). We 
hope that these predictions will act as a spur to further 
experimental work. 
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